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Abstract 

We use the symmetry of c = 1 quantum gravity to compute matrix model special state 
correlation functions. The results are compared, and found to agree, with expectations from the 
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1 Introduction 



The last few years have seen tremendous developements in the understanding of two dimen- 
sional quantum gravity and therefore non critical string theory. The first success came from 
field theory, |l|, where the case with matter of central charge c < 1 coupled to quantum 
gravity was solved. Later these theories were also solved by the powerful matrix models, 
both for c < 1 0, and for c = 1 p|. They allow for exact, non perturbative solutions where 
one sums over all genus. This is an important improvement over field theory where higher 
genus contributions are extremely difficult to calculate. Most results have therefore been 
limited to the sphere or in some cases the torus. 

These models are important for two reasons. First they give solutions of two dimensional 
quantum gravity and can serve as toy models for higher dimensional theories. Indeed, the 
matrix model provides directly the exact Wheeler de Witt equation summed over all genus, 
HJ]. Clearly an important object for further study. Second, and perhaps even more inter- 
esting, they describe non critical string theories. The most physical example is clearly the 
case c = 1. It is now commonplace to identify the Liouville mode as an extra dimension, ||, 
and thereby obtaining a theory of strings moving in a two dimensional target space. Naively 
one would expect this theory to be very simple, just a single massless scalar particle, usually 
referred to as the tachyon. This can be argued by choosing lightcone gauge. Fortunately, 
this is not the whole story. Instead there are remnants of the massive excited string modes 
present for certain discrete values of the momentum, ||, ||, [7], ||. They are usually called 
special or discrete states. Clearly it is very interesting to study these extra states if one 
wants to learn about truly stringy phenomena. 

The special states have been the object for several recent studies. Both using Liouville 
theory and matrix models. An important discovery has been a huge set of symmetries. These 
symmetries obey a algebra, which can be thought of as a generalization of a Virasoro 
algebra. This has been shown in the Liouville theory with two different methods. In || the 
Woo was found by explicitly calculating the operator product expansion of the special states. 
An important tool for doing so is the usual SU(2) symmetry known for a long time. In 



10] the symmetry was instead found from the construction of the ground ring. This gives a 



representation of the current algebra with the currents acting on a set of ghost number zero, 



spin zero fields. These fields can be shown to be primary, |TlJ, and correspond to special 
states in addition to the standard ones. The meaning of these and other special states of non 
standard ghost numbers has as yet not been fully clarified. From the matrix model point of 



view, the emergence of the symmetry has been more gradual ]T0, 12, 13, 14 1. In [ 10 1 it was 



however clearly realized that the simply is generated by the matrix eigenvalue and its 
conjugated momentum through the Poisson bracketts. 
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In this work we will use the to study the special operator correlation functions. 
Notations and conventions will be much the same as in a previous paper, ||. We will be 
able to obtain many of the results in || using the simplifications the symmetry provides. 

In section 2 we make some initial comments on Ward identities and symmetries relating 
to the results obtained in ||. We will also make some comments on possible generalizations 
to non harmonic matrix model potentials based on a generalized Wheeler de Witt equation. 
In section 3 we calculate some matrix model special state correlation functions using the 
symmetry. We also make some comments on how to identify the counterparts of the Liouville 
model special states. Section 4 gives some illustrations of the structure one encounters at 
higher genus. Finally, section 5 is devoted to a comparison with Liouville theory. Although 
the success of the matrix model and its agreement with Liouville theory hardly is in doubt, 
it is important to make the connection as explicit as possible. In particular, it is so far 
not clear how to explicitly extract the space time structure from the matrix model. One 
would like to be able to study nontrivial space times like the recently discovered black hole 



solution We will not be able to adress this question here, but we will be able to make 
a comparison between our matrix model results and some Liouville theory expectations. In 
the cases which we will examine we will find perfect agreement. 



2 Ward Identities and Symmetries 

Let us now consider the matrix model and its special states and operators. The matrix 
model represents the string theory Riemann surfaces by Feynman diagrams of interacting 
matrix variables which triangulate the surfaces. In the uncompactified case, or at least for 
large enough radius, we can simply integrate out most of the degrees of freedom. The only 
remaining will be the matrix eigenvalues. They will then behave as non interacting fermions 
in the matrix model potential. Obvious candidates for special state correlation functions 
are then correlation functions of powers of the matrix model eigenvalues. Such objects were 
studied in || and the expected poles for discrete momenta were found. In this section we 
will review and extend some of the results of || for Ward identities. These Ward identities 
can be used to recursively obtain the correlation functions. 

The recent developements revealing the symmetries indicate however that this is 
not the whole story. One should also consider correlation functions involving powers of the 
conjugate momentum. In this section we will show the existence of this symmetry which 
greatly will simplify the subsequent calculations. 

Ward identities for correlation functions are in general obtained by changes of variables 
in the path integral. Examples of such Ward identities were obtained in [H] from simple 
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coordinate changes in the matrix eigenvalue. They can be thought of as generated by com- 
mutators, or classically i.e. on the sphere, Poisson bracketts with pX m . They obviously obey 
a Virasoro algebra which is part of a algebra generated by all monomials p n X m . We 
may also introduce time dependence and consider generators with certain momenta q, i.e. 
p n \ m e iqt . Let us as an example make a pX k e iqt variation of the one puncture function. The 
two puncture function is schematically given by 



DC 







<PP>=lm dT dXG(X,X;T) (1) 



where the calculation is done at the Fermi surface. We will shift its energy to zero, hence 
putting the Fermi energy as a constant term in the potential. G is the path integral given 
by, in Euclidean time, 

G(Ai,A 2 ) = [ X2 [dpdX]e-^ dt ^ p2+u ^ (2) 



where U(X) = ^ p t p X p is the potential. The variation of the partition function would have 
involved a sum over all states in the Fermi sea up to the Fermi level. By inserting a puncture, 
i.e. taking a derivative with respect to the Fermi energy, we restrict ourselves to the Fermi 
surface. Next we perform the variation of < P >. The measure, as given by (0), is invariant 
under the change of variables. (Clearly we are not supposed to differentiate with respect to t 
when changing variables in the measure.) The change in the action give rise to the following 
identity among two point functions: 

r dU 
< ( J dt(iqX k p + kX k ~y + —X k )e iqt )T >= (3) 

The puncture is now a tachyon, T, carrying away the momentum. The piece with a single p 
is evaluated by integrating over p in the path integral obtaining a A which then is partially 
integrated. We then switch to a Hamiltonian formulation, remembering that we should use 
Weyl ordering. We finally obtain: 

- -k(k -l)(k-2)< O fe _ 3 T > q , g -i + ^2(2k + p)t p < Op+t^T > q>g (4) 
1 v 

where q indicates momentum and g genus. We have introduced the notation Ok for X k . The 
different genus for the first term is due to a p, X commutator which arises when we want to 
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evaluate the p 2 against the wave function of the Fermi surface. This gives an H, which is the 
same as the genus coupling constant. Following |4j] we may define the loop operator given 
by 

w{l) = e lx (5) 

It corresponds to cutting out a hole in the surface with a boundary of length I. The reason 
is as follows. If we insert a power n of the original matrix eigenvalue m on the surface this 
creates a little hole, the length of the boundary being proportional to n (the number of legs) 
and the lattice spacing a. For fixed n the length clearly shrinks to zero in the double scaling 
limit. To get a finite length we must also take n to infinity. Introducing A as m expanded 
around the top of the potential, we find in the double scaling limit 

m n ~ (1 + a\) l/a -> e lx (6) 

We may then Fourier transform to obtain a differential equation in the loop length. We get 

E + + to12 + q2] < w{l)T >= (7) 

where the third term are of order one higher in the string coupling and do not contribute on 
the sphere, to is the Fermi energy \x with the appropriate number of /3's absorbed. In the case 
of the usual harmonic oscillator potential, where to = fin, the resulting equation is in fact 
the Wheeler de Witt equation obtained in |4]] with a more indirect matrix model method. 
On the sphere this is one of the most striking verifications of the equivalence of the Liouville 
and matrix models. At zero momentum the Wheeler de Witt equation is just the Fourier 
transform of the Gelfand-Diiki equation for the resolvent for the Schrodinger operator. This 
was the way in which the zero momentum version of (f4|) was derived in . If we want to be 
careful, see section 5, we need to rescale A by y/—t 2 to get a dimensionless I. This is needed 
for the exact correspondence between the above result and the mini superspace canonical 
quantization of Liouville theory. Recall that, [0], t 2 = — J^y. From (0) one might try to draw 
some conclusions about the Liouville theory correspondence to the more general potentials 
above. Clearly the last term, which corresponds to the matter piece, does not change while we 
change potential. Instead it is the piece which would be expected to arise from a canonically 
quantized kinetical term for the Liouville mode which gets modified. Hence one is lead to 
the conclusion that these more general models (however with p independent potentials) may 
correspond not to modifications of the matter theory but rather to different theories for the 
Liouville part. This is also consistent with the point of view for which this paper will argue, 
that the special states must be represented using both A's and p's, not just the A's. 
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In principle all correlation functions may be calculated with the help of Ward identities 
derived in this way. However, it is more convenient to make use of the large set of symmetries 
in the theory. As shown by Witten in jlTJ we may change basis to (p — X) n (p + A) m and, 
for certain time dependence, obtain transformations which leave the action invariant. These 
transformations are generated by, in Minkowsky time, 

W r ' s = (p + X) r (p- X) s e (r - s)t (8) 

Again we get the algebra 

{W ri ' Sl , W r2 ' S2 } = (ns 2 - r2Sl )M/ ri+r2 - 1 ' sl+S2 - 1 (9) 

generated, classically, by the Poisson bracketts. For a general momentum q in (H) we find 
when acting on the Minkowsky action S = J (pX — \{p 2 — A 2 )) 

{W r ' s , S} = (r-s- q)W r ' s (10) 

Hence a symmetry for appropriate discrete values of imaginary momentum. This is equivalent 
to saying that W = W(p, A, t) is a solution of 

dW dW 

_ = _ + {JW = (11) 

Expressed in terms of the initial conditions, po and Ao, we have W = W(po,Xo), i.e. any 
time independent function of the initial conditions. The generators (El) are then simply 
obtained through evolution in time. Hence the transformations can be understood as time 
independent canonical transformations of the initial conditions. 

A more convenient way of labeling the generators is through their SU(2) quantum num- 
bers J = (r + s)/2 and m = (r — s)/2. With the definition 

Wj, m ={p + X) J+m {p - X) J ~ m e 2mt (12) 

one gets 

{Wj umi ,Wj 2>m2 } = 2(mi J 2 - m 2 J 1 )W Jl 

+J2-l,mi+m 2 

(13) 

In Euclidean time, which is what we will be using, one should take p — *> ip and t — > it in 



121) . There is also an extra i in the structure constant of (p~3|) and the eigenvalue of (|T0|) . 
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3 Matrix Model Correlation Functions 



In this section we will calculate correlation functions involving the operators W as defined 
above. The Woo symmetries will help us organize the Ward identities. As an example, let 
us start with the two point function. Using 

1 00 1 

<PP >=-Im£ (14) 
and simple perturbation theory we get: 

< W 1 W 2 P >= -Im 2^ 2^ (15) 

71 timeord n,k E ™ + £ ° ^1 + E k ~ E n 

Since the W's are of the form (§), continued to Eucledian time, they are simply raising or 
lowering operators in the inverted harmonic oscillator. This means that only a few of the 
matrix elements are nonzero. Since W raises by 2m = r — s we get 

< WlWa P >= l ImE <H[^ 2 ]|n> 1 

We have reduced the two point function to a one point function using the commutation 
relations. If we restrict ourselves to the sphere, use the algebra given by @ (with an extra 
i in the structure constant for Eucledian time) and directly calculate the one point function 

we get 

< Wl W 2 >= 2(^-^) 1^1 , , ogA , 2 , 1+ , 2 -, (17) 

2mi - pi 7r Ji + J 2 

or equivalently 

< WiW 2 >= — — -/i Jl+Ja I log// I 2 Jl+Ja (18) 

2m! — pi n 

The simplest way to obtain the one point function on the sphere is to use the classical 
Fermi liquid picture introduced by Polchinski JHJ. We simply need to do the phase space 
integral: 

<W >= J dpd\W(p, A) (19) 

This was also discussed in [|HJ . To make everything well defined we need however to introduce 
an extra puncture, i.e. take a derivative with respect to the cosmological constant. Doing 
that the integral over the whole Fermi sea becomes just an integral over the Fermi surface: 

< WP >= <fw(p,X) (20) 
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The integral is to be performed along a hyperbola p ~ A = /t. Although the answer is really 
infinite, we know that the piece nonanalytic in \i is - | log/i | for < PP >, i.e. when 
we just integrate 1. The other zero momentum operators simply involve integrations of 

2 \ 2 

^' 2 - = — ^ n = /i™, again constants along the Fermi surface hyperbola, so we get 
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< W n ' n P >= -fi n I log/i I 2 n (21) 



7T 



and from this ([17]) follows. Our conventions are such that a' = 1. Another approach, which 
is convenient when calculating correlation functions of nonzero momentum, is to continue 
to the upside down oscillator where the Fermi surface is a circle. In that case, however, we 
need to remember to put the Liouville volume | log fi | in by hand. Parenthetically we may 
note how a general correlation function may be obtained in this way. For instance the two 
point function is obtained by perturbing the hamiltonian and hence the Fermi surface by 
one of the operators. If we integrate the other operator against the change in Fermi surface 
we get the correlation. 

Let us now consider the more complicated case of a three point function. Again pertur- 
bation theory gives us 

1 x , ^ <n\Wi\m X m\W 2 \k X k \ W 3 \ n > 1 1 

— Im > y 

71 Uncord n^fc E n + *0 l Pl + E m - E n tp 3 + E n - E k 

(22) 

Let us make the sum over time orderings more explicit. We find 

1 <n\(W 1 W 2 W 3 + W s W 2 W 1 )\n> 1 

* „ K+k i(pi - 2m 1 )i(p 3 - 2m 3 ) ' ( ' 

This may after some straightforward manipulations be rewritten as 

< Wl w 2 w 3 >= UgkElM> + lPMg^Jl> m 

(pi - 2m 1 )(p 3 + 2mi + 2m 2 ) (P2 - 2m 2 )(p3 + 2m x + 2m 2 ) 



For the sphere we now use the algebra (|13| ) and the explicitly calculated one point function 

to get 

_ (2mi - pi)m 2 m 3 J 1 + (2m 2 - p 2 )mim 3 J 2 + (2m 3 - p 3 )m 1 m 2 J3 

< kV \W2Vv 3 ^* — 7 77 77 : 

(2m 1 - pi) (2m 2 - p 2 ) (2m 3 - p 3 ) 



x-/x^ - 7 "" 1 I log// I 2^"=i J ™ (25) 

7T 



The general higher point function can be obtained recursively from the three point by 
use of (HD and (|10|). To get the N point function with an additional operator Wn we vary 
the N — 1 point function with Wn knowing that the total variation is zero. The variation 
consists of two terms. One from varying the action as given by fllPp and a sum of terms from 
varying the other operators as given by (|13|). Each of the terms in this sum is obtained by 
shifting pi — > Pi+PN, mi — > m>i + and Jj — > Jj + Jn+i — 1- We also need to multiply with 
the Clebsch-Gordan coefficient 2(^171^ — Jn^i)- It is an easy exercise to check that (25) is 
obtained by applying this procedure to ([[7]) or ([18]). There is one subtlety in the variational 
procedure which should be noted. The insertion of an operator W in the pathintegral does not 
involve only the operator itself, but also a delta function for its position in eigenvalue space. 
In general the delta function will also contribute to the variation. Luckily its contribution 
will be zero by invariance properties for the Woo generators. 

Using this method one can write down several different recursion relations. One simple 
example is: 

< T Jt jWj 1>mi n >= 2J( /;~ mi) < wj +Jl ^j +mi n T JitJt > ^ 

i=2 ZJ P i=2 

We will come back to this relation later, when we compare with the Liouville model results. 

Rather than considering these general expressions, let us look at a couple of important 
examples where the form of the general N point function is particularly simple. 

The first example is the N point function of special tachyons. It is given by 



Jl 2 T\ N i Jr, 1 ^71 1 

< n r « >= ^ r T - c ro^- 1 " i kg** 1 2*:- j - (27) 

n=l 



n^ 2 (2Jn-Pn)vT^- 3 



The quantum numbers have been chosen as m n = J n for n > 1 and m\ = —J\. This is just 
the pole part of the general tachyon correlation function as computed both in the matrix 
model |18| and in the Liouville theory |18, 1^], up to a factorized normalization factor. The 



proof is by varying the three point. We can not just vary the three point tachyon correlation 
function, since some of the J's are really m's in disguise and J and m vary differently. Instead 
we start with the general three point and make an arbitrary number of tachyon variations. A 
simplification is that we at each step only have to vary the single negative chirality tachyon. 
It is only from there were we will get a nonzero Clebsch-Gordan coefficient. Following the 
prescription above, performing iV — 3 variations we get 
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N N N 

< II T n >= [(2mi -pi + E( 2m * - Pi))m 2 m 3 ( J x + ]T J % - iV + 3) 

n=l i=4 i=4 

AT AT 

+ (2m 2 -p 2 )(m 1 + ^mi)m 3 J 2 + (2m 3 -p 3 )(mi + E m *) m 2^] 

x n£l^ + gfe±l^j - ( J i + g|i±l 4 - g + 0^) j ^ i j i 2 E1 1 ^(28) 

(2m 1 -p 1 + E£ 4 (2mi-p l ))nn= 2 (2m n -p n ) v 1 &P 1 v ; 

The product in the denominator is the product of all the Clebsch-Gordan coefficients of the 
variations. Note that each get shifted by the successive variations. By the use of momentum 
conservation and evaluating the m's as J's, the formula (p7|) is proved. This derivation shows 
how the combinatorical factor from the \x derivatives is a consequence of the symmetry. 

If we want to consider the zero momentum operators, we have to be careful. The Clebsch- 
Gordan coefficients are zero in this case but these zeroes cancel precisely the momentum poles 
and leave a finite result. Also, we need to consider both signs of the m quantum number 
when we take the m — > limit. This gives a necessary extra factor of two. We get 



N N i JN-3 N N 

< II W n >= E Jn-T^f^^' 1 I logj* I 2En =1 ^n (29 ) 

71=1 71=1 71 a ^ 



We will use induction for the proof. We find 



AT+1 ^pN+l j _ -1 

II li v >= X)2(Jjv +1 m fc - •//,/// v., ; 1 ' " 



„=1 k 2mAr + i - p N+ i 

x ~ dfI N-3 ^ n=1 I lo SA* I 2^™= 1 n + (wat+i -> -miv+i) (30) 

If we then put Pat+i = and use that the sum of all m's must be zero the result follows. 
This can also be checked by an explicit phase space calculation. 

Given these expressions we may check the correlation functions calculated in || and 
independently in [|]]. These were correlation functions of pure powers of the matrix eigenvalue 
without any momentum powers. In terms of the Ws they are given by 

n = E ( I ) (-1)^-^ (31) 
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From this it follows that the two point function is given by 



< O n O m > q = ^ EE [ k )[ i ) (-±) k+l < Wr k < k W- 1 ' 1 > (32) 
Using (^) and some simple algebra we find 

-J—ljM\ ]o zu\Y( n )( m ) 4(1 ~ k? (S3) 

2 I(n +m) ^ 1 l0g/i 1 ^ \k ) \ 2±2 - k ) 4(f - kf - {66) 

In precise agreement with || recalling our convention a' = 1. We can now understand 
why the O operators gave correlation functions with sets of poles and were, depending 
on momentum, capable of exciting several special states |J. They were, in fact, linear 
combinations of all special operators of a given gravitational dimension i.e. spin J. The 
above construction with the generators @ of the Woo disentangles the correlation functions. 
This means that the matrix model operators to be identified with the Liouville model special 
states are those defined in ([12]). 

Finally let us consider the meaning of the momentum poles. As emphasized in |23] we 
should not treat the poles in (25) and the | log/x | asymmetrically since the source of the 
| log/i | is also a momentum pole. In fact, all the poles should be thought of as cut of by 
| log fi | . A general iV point function (without zero momentum operators) would then have 
| log/i 1^. This proliferation of logaritms was also noted in W9 . 



4 Higher Genus 

We have so far basically just treated the sphere, which means, in the matrix model, that we 
have been working at the classical level. The has been generated by Poisson bracketts. 
Nothing can however stop us from considering the full quantum theory, i.e. all genus. It is 
just a matter of algebra to compute for instance the two or three point functions using (|16D 
and (0) respectively. 

More interestingly, the algebra changes at the quantum level. There is a deformation with 
h, the genus coupling, as parameter when we use commutators instead of Poisson bracketts. 
If we define our Ws using Weyl ordering, which is natural from the pathintegral point of 
view, the algebra may be conveniently represented using the Moyal bracketts |20| 

{WuW2] ^l s ^A. r A.A. )WlW2 ,34) 
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which is a deformation of the usual Poisson brackett. From this we might conclude that also 
the Liouville theory operator product expansions should receive higher genus corrections. 
Presumably from handles getting caught inside the contour integrals defining the operator 
product expansions. 

Another way to exhibit the quantum deformation is through a generalized loop operator 



lq| , instead of (g) we introduce 

w(k,l) = e kp+lx (35) 



where p and A are the conjugate variables. It is a very old result, that these operators 
obey the algebra 

[w(ki, h),w(k 2 , h)} = ~ sin %{k\l 2 - hki)w(ki + k 2 , l\ + h) (36) 



with % — > giving back a W^. Interestingly it can be shown |2T| that (|3q) is a representation 
of SU(N) with H = 1/N. This is reminiscent of the original unitary symmetry of the matrix 
model. 

Let us give an explicit example of a two point function to all genus. We choose the 
correlator between spin J = 3/2, m = 1/2 and J = 3/2, m = —1/2. To do that we need to 
calculate < W^flP >■ This is easy. We have 

<iy 20 P>=-ImV {2E ^ 2 + l (37) 

The extra term +1 comes from Weyl ordering. If we keep only terms nonanalytic in to this 
reduces to 

< W 2fi P >= {Atl + 1) < PP > (38) 
To evaluate our two point we use the Moyal brackett to calculate 

[^3/2,1/2, W3A-1/2] = 6^ 2 ,o - 4h 2 W 0fi (39) 
(|I6D and (|38| ) then finally give 

< W 3/2A/2 W 3/2 ,„ 1/2 P >= (48tg + 4) < PP > i^- (40) 
The same procedure may be used to calculate arbitrary correlation functions. 
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5 Comparison with Liouville Theory 



We would like to understand the structure from the Liouville theory point of view. As 
shown in |10| and by more direct methods in ||, we indeed have the same algebraic structure 



present. Therefore one would expect the comparison of the Liouville and the matrix model 
to be straightforward. As we will see, the situation is more subtle. Let us first consider a 
very special case which also give us the opportunity to clarify some important points. 

There are some very simple examples of correlation functions easily computable just 
using Liouville notions and no matrix model techniques. These are correlation functions 
involving the dilaton. We will in fact be able to obtain some results to all genus simply from 
dimensional arguments. Consider the Liouville partition function (or space time free energy) 

E(A) = lim - [ DX D(f ) e-f^ dx8x+d ^ + ^ +Aea ^ (41) 

R^co R J 

where t 2 = — 7^7, Q = 2y/2, a = —\/2 and R is the radius of the target space for the matter 
field X. A is the world sheet cosmological constant, dimensionless from the point of view of 
space time. The only dimensionful quantities are R and a'. In the noncompact case we have 
infact only a' at our disposal. From dimensional grounds and KPZ scaling we must have 

E(A) g ~ (-t 2 )V2A 2 ( 1 -fl) (42) 

at genus g. E(A) is the generator of connected amplitudes (in space time). Let us do a 
Legendre transform to obtain a generating functional for 1PI amplitudes with respect to the 
puncture, i.e. the zero momentum tachyon. This means taking away any pinches. We have 

E(A) — Afj, — TQm) (43) 

with 

A*=|f = (-t a ) 1/2 (A + ...) (44) 
/j, has the dimensions of energy. Hence 

r(ji) ~ (-fc^-V -2 ' (45) 

In the 1PI generator T we should of course regard \x as independent of t 2 . Since t 2 derivatives 
should generate dilaton insertions we find the following 1PI amplitude relation 

1 n 

< 2 ...0 2 > g = — — U(2g + 1- 2p) <> g (46) 
\ M V P =l 
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which is identical to what was obtained in || using the matrix model recursion relations 
(generalizations of the zero momentum Wheeler de Witt equation). As noted in || the 
dilaton one point function involves a factor 2g — 1 rather than the expected 2g — 2. From 
above it is clear that this discrepancy is simply due to including the overall (— 1 2 ) 1 ^ 2 in ([42]). 
There is a further subtlety in how the dilaton is defined. As we have seen a pure A 2 is not 
what we would expect to identify with the dilaton. Instead we should have p 2 — X 2 if we 
keep the algebraic structure in mind. These are in general different in correlation functions. 
We we will return to this shortly. 

It is important to realize the difference between connected and 1PI amplitudes. In the 
matrix model 1PI amplitudes are the natural objects, in Liouville theory it is more common 
to treat the connected ones. Often the distinction is not very clearly made. Indeed if we 
consider generic nonzero momentum the difference is very easy to deal with. It amounts 



to a renormalization of the cosmological constant. We can simply replace A by /i, [23]. At 



zero momentum we must be much more careful. In this case we may have internal puncture 
propagators, i.e. pinches. Some examples of this were obtained in |J. This will turn out to 
be important later on. 

Already at this point we may find traces of the structure. In fact, the seemingly 
innocent representation of the puncture and the dilaton as \i and t 2 derivatives respectively 
is a reflection of the W^. Let us give a formal argument for this. First the puncture. Write 
the SU(2) quantum numbers of the puncture as J and m which both will be taken to zero. 
Choose one of the operators in (25) to be a puncture. We get 

< W X W 2 P >= Ui + J 2 ) — — -/i Jl+j2_1 I log// I 2 Jl+Ja (47) 

2mi — pin 



which by comparing with the two point function (|r7|) shows how the puncture is represented 
as a n derivative. The case of the dilaton is equally simple. Proceeding as above we find 

< W,W 2 D >= [(J, + J 2 )-^- + - 2mj d V 1+Ja I log// I 2 J ^ +1 (48) 

2m\ — pi (2m 1 — pi) 1 7i 

If we introduce explicit t 2 s in the two point function we can write it as 

< W X W 2 >= ^ 7T77-(t — ^tt7^) Ji+J2 I log// | 2 Jl+j2 (49) 

2m 1 -p 1 /(-2t 2 y/ 2 7f\-2t 2 ) 1 / 2j 1 ^' K J 

Taking a t^ 2 derivative we indeed reproduce fl48|) . We must now return to the issue of how 
precisely the dilaton is defined. Recall the original matrix model action: 

(3 J dt[p\ - l -p 2 - t 2 \ 2 } (50) 
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with f3 dimensionless, ti having the dimension of energy squared and p 2 and A 2 the dimensions 
of energy and one over energy respectively. To obtain (|49|) as a generating functional for 
dilaton insertions with the above definition of the dilaton we should rescale A and p to make 
them dimensionless. We find 

(3 J dt\p\ - -^=(-t 2 ) 1/2 (p 2 - A 2 )] (51) 

Hence the matrix model dilaton should be represented by (— t2) 1 ^ 2 derivatives. This is the 
rescaling eluded to in section 2 in the context of the Wheeler de Witt equation. 

To obtain the general special operator correlation function in the Liouville theory one 
would like to use the group theoretic information provided by the or, for given spin J, the 
SU{2) symmetry. The states in the Liouville theory are given by combinations W(z)W(z) 
of the Liouville theory version of the special states W(z). Given this it is tempting to believe 
that we have a representation of a x Woo symmetry (left times right). This is however 
in general not correct. In the uncompactified case the left and the right moving states must 
be the same. The symmetry group is broken down to just the diagonal subgroup. This is 
achieved in two steps. First the gravitational dressing must be the same for left and right, 
otherwise we would be unable to screen using the cosmological constant which treats left and 
right in the same way. This means that we always must have the same spin J for left and 
right. We get a reduction to the diagonal of the piece transverse to SU(2) x SU{2). This 
is true even for the compactified case. If we in addition are considering the uncompactified 
case, the left and right moving momenta must be the same and hence the m quantum 
numbers. Consequently we just have a representation of the diagonal Woo- This is in precise 
agreement with the matrix model, where we indeed only see one Woo- There is however 
an apparent paradox here. If we would use the free field contractions in computing the 
correlation functions the results would seem to disagree since from this point of view left 
and right are still independent. We will return to this important point further on, and 
discover that there in fact seems to be no contradiction. 

The symmetry may then be used to determine all correlation functions given the special 
tachyon correlation functions which may be computed using other means. The reason is 
that all J and m dependence of any correlation function is given by some combination of 
Clebsch-Gordan coefficients. For given J's we need the Clebsch-Gordan coefficients of SU(2), 
the 3j symbols, to get the m dependence. In fact we have already seen the agreement for 
the tachyon correlation functions and if we believe that the group theoretic structure is 
the same in the matrix model and in the Liouville theory, we know that the expressions 
obtained in the matrix model must agree with Liouville theory. To be more explicit let us 
however look at an example, the three point function, to see how the invariance properties 
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determine the correlation functions. The three point function is obtained by considering 
coupling ( Ji, mi), ( J 2 , m 2 ) and ( J 3 , m 3 ) (with m 1 +m 2 + m 3 = 0) to ( J\ + J 2 + J 3 - 2, 0). A 
complication is that there are in general several different channels to sum over. This is true 
already for the three point function. The reason is that we really should think of the three 
point function as a four point function. The fourth leg carries the excess Liouville momentum, 
i.e. J quantum number, into the vacuum. This is a consequence of the non conservation of 
Liouville momentum. Let us use the tachyon three point function for normalization. It is 
given by 

< T{T 2 T 3 >= ^ I log/. I 2*+*+ J »+ 1 (52) 

(2m 2 -p 2 ){2m 3 - p 3 ) n 

where we have kept the normalization choice of fl2"7|). Tachyons 2 and 3 are of positive chirality 
while tachyon 1 has negative chirality. There are two possible channels corresponding to 
either p 2 = m 2 or p 3 = m 3 , i.e. 1 and 2 coming together or 1 and 3 coming together. The 
group theoretic factor in each case is simply proportional to a product of 3j symbols. One 
for each vertex. For the 1-2 channel: 



JiJ 2 J \ / JJ 3 J' 
m\m 2 m 3 J \ —m 3 m 3 



(53) 



where J = J\ + J 2 — 1 and J' = J\ + J 2 + J 3 — 2. Just retaining the m dependence and 
adding the two channels we find 

(2m 2 - p 2 ){Jim 2 - J 2 nii)m 3 + (2m 3 - p 3 ){Jiin 3 - J 3 m 1 )m 2 = 
-(2mi - pi)Jim 2 m 3 - (2m 2 - p 2 )J 2 mim 3 - (2m 3 - p 3 )J 3 m 1 m 2 (54) 



which agrees with (25) after using ([52]) to fix the normalization and J dependence. This 
should come as no surprise since the calculations are almost identical. 

Another convenient way to obtain more general correlation functions is through factor- 
ization. This is really already implicit in our previous calculations. In fact, if we look at 
( P5| ) we see the complete factorization of the tachyon correlation function into a product 
of three point functions, each given by a 3j symbol, times a single zero momentum one 
point function. This last piece represents the extra leg in any correlation function which 
absorbs excess Liouville momenta. One may note that these three point functions in fact 



involve states of the wrong dressing. This was also pointed out in p4fl . Strictly speaking 
the expression in (|28| ) is just for one channel, the one where 1 fuses with 2 then with 3 etc. 
All channels give however identical contributions and can not be distinguished. Clearly the 
tachyon correlation function is consistent with the single factorization result. 
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As has been remarked, this seems to be in contradiction with what to expect from the 
naiive free field calculations in Liouville theory. From such a calculation you would expect to 
get a different result, all Clebsch-Gordan coefficients squared, one from the left and one from 
the right. We will however show that the results in the end turn out to be consistent. Let us 
begin by considering the tachyon correlation function as computed in [19|. As we have seen 
the result is in complete agreement with the matrix model results. On the other hand we 
have seen how the matrix model organizes its correlation functions using a single W^. Let 
us consider the Liouville calculation more carefully. The result of |1| is obtained through 
arguments of analyticity and symmetry. In particular the by now well known factorized 
product of gamma functions is found with a certain unknown coefficient independent 
of the particular momenta. This coefficient is then determined by sending all the momenta, 
except three, to zero. This reduces the expression to a three point function with N — 3 
extra punctures. Since the three point function is possible to evaluate directly, the general 
result follows. The extra N — 3 punctures is simply represented as /i derivatives. This is 
the Liouville derivation of the expression (|2~T|). The important point is that the use of a /i 
derivative for inserting a puncture is a consequence of having just one Wool This means that 
the calculation in [19|] automatically incorporates this feature. 

For the more general case with nontachyonic special states, we return to the recursion 
relation (2G). Let us redefine the fields according to 

2J ~ 

Wj, m = Wj, m (55) 

2m — p 

Then the recursion relation takes the form 

< f JtJ w Jumi ft f JuJi >= (Ji ~ mi) y +Ji ~ 1} < www+ mi n Tj itJt > (56) 

i=2 J l i=2 



In p5[ these very same recursion relations were obtained in the case J = m = 1/2 using 
Liouville methods. The coefficient in front of the right hand side were shown to be of the 
form (2Ji — 1)C 2 , where C stands for the appropriate Clebsch-Gordan coefficient. The first 
factor comes from comparing with the purely tachyonic case where the answer is obtained 
from a simple Veneziano like integral. To see the agreement one uses the Clebsch-Gordan 
coefficients of the special operator algebra as obtained in [[J . 

= w^) < j ™ " (57) 

where 

A(J,m) = -i[(2J)!(J + m)!(J-m)!] 1/2 (58) 
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At J = m = 1/2 one finds C 2 = which then leads to This is an important check 

on the equivalence between the Liouville and matrix model approaches. An everywhere 
present difficulty in these comparisons is, however, the fact that we are really sitting right 
on the momentum poles. Clearly one needs to carefully regularize all expressions. 

Let us give some further illustrations in the case of puncture and dilaton insertions. 
We begin with the puncture. Starting with a general correlation function and inserting a 
puncture does not change the Veneziano like integral which has to be calculated. When we 
insert a puncture we also must remove one of the screening insertions. The only thing which 
changes is the zero mode part of the calculation. We recall the result 

4>e m ^ Ae ~* ~ T(-m)A m (59) 

If we start with T(— m)/i m we end up with T(—m + l)/i m_1 = — mT(— m)/i m_1 when we 
remove a screening insertion. Comparing with p7| ) this shows the origin of the related 
factor in the tachyon correlation function. It is a consequence of the changing number 
of puncture screening operators needed. For the tachyons the issue of connected or 1PI 
amplitudes is trivial for our case with just one tachyon of differing chirality. There can't be 
any internal punctures just from kinematics. This is no longer the case when we turn to the 
dilaton. The crucial point is that the dilaton can be represented as a t 2 derivative. Usually 
this is precisely equivalent to using the ordinary free field contractions giving Veneziano like 
correlation functions. Inserting a dilaton in some tachyon correlation function means taking 
derivatives with respect to t 2 (i-e. 1/a'). For dimensional reasons all tachyon momenta are 
accompanied by an t 2 . Without explicit t 2 s one could write k-^ for the dilaton. For a dilaton 
insertion in a nonzero momentum correlation function the dominating pole contribution 
comes from letting the t 2 derivative act directly on the poles. All other terms are clearly 
less singular. This gives the second term in (fy|). At zero momentum we must also consider 
the dependence from the t 2 s which go together with the /z's. The latter is a consequence 
of dealing with 1PI rather than connected correlation functions. One could in fact obtain 
the result by considering the explicit combinations of 1PI amplitudes into connected ones. 
To be more precise, if we want to obtain the 1PI amplitude from the connected amplitude, 
we must amputate external puncture legs but also subtract of the diagrams with internal 
puncture propagators. In particular we need to subtract a diagram where a puncture goes 
off and converts into a dilaton. This diagram therefore involves a puncture insertion and 
gives a contribution corresponding to a fi derivative. This is then simulated by an explicit 
t 2 accompaning the /i's to assure the proper subtraction. This corresponds to the first term 



in (f48l) , which only becomes relevant compared to the second term at zero momentum. 



Otherwise we will get one power less of | \ogfi |'s. We need not restrict ourselves to a 
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dilaton among special tachyons, the same reasoning works for a dilaton inserted in a general 
correlation function. From these two examples we can conclude that the 1PI nature of the 
correlation functions is very important in the case of zero momentum. 

6 Conclusions 

We have investigated the structure of special operator correlation functions in c = 1 quantum 
gravity. Due to the presence of a Woo symmetry the calculations become very simple. We 
have also investigated the connection between the Liouville and the matrix model, indicating 
the agreement for the correlation functions. 

An important point is the existence in the uncompactified Liouville as well as matrix 
model formulation of c=l of just one W^. From the Liouville point of view this is somewhat 
obscure since the operator product expansion and its Clebsch-Gordan coefficients seem to 
give a structure corresponding to two W^s. Fortunately the final outcome of the explicit 
calculations are identical. Further work is however needed to establish the full equivalence. 

A part of the problem is that many of the calculations are so ill defined. The reason is 
that we are sitting right on the discrete momentum poles. Especially in the Liouville theory 
this is a big technical problem. Often we must rely on guesswork concerning ill defined 
analytical continuations. It is very doubtful if many of the results would have been obtained 
correctly without knowing the answers in advance, given by the much more powerful matrix 
model. 

Important issues for future research are to investigate multicritical points of matrix mod- 
els with generalized potentials. It is natural to consider even non quadratic dependence on 
the momentum. Such perturbations could arise from adding the special states as we have 
seen. Another important point is to identify the 'wrongly' dressed special states in the ma- 
trix model. Such states have negative gravitational dimensions and are important in the 
context of the two dimensional black hole. 
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